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Abstract. The generalised longitudinal susceptibility xil^'^) affords a sensitive 
measure of the spatial and temporal correlations of magnetic monopoles in spin ice. 
Starting with the monopolc model, a mean field expression for x(*il) ^) derived as 
well as expressions for the mean square longitudinal field and induction at a point. 
Monopole motion is shown to be strongly correlated, and both spatial and temporal 
correlations are controlled by the dimensionless monopole density x which defines the 
ratio of the magnetization relaxation rate and the monopole hop rate. Thermal effects 
and spin lattice relaxation are also considered. The derived equations are applicable 
in the temperature range where the Wien effect for magnetic monopoles is negligible. 
They are discussed in the context of existing theories of spin ice and the following 
experimental techniques: dc and ac-magnetization, neutron scattering, neutron spin 
echo, and longitudinal and transverse field /uSR. The monopole theory is found to unify 
diverse experimental results, but several discrepancies between theory and experiment 
are identified. One of these, concerning the neutron scattering line shape, is explained 
by means of a phenomenological modification to the theory. 



Generalised Longitudinal Susceptibility for Magnetic Monopoles in Spin Ice 



2 



1. Introduction 

Following the paper of Castelnovo, Moessner and Sondlii [1] on emergent magnetic 
monopoles, there has been renewed interest in the properties of spin ice O |3l HJ |5]. 
Magnetic monopole currents were first envisaged by Ryzhkin |£j while Jaubert and 
Holdsworth [3, E] studied the closely related problem of magnetic relaxation in spin ice 
by means of numerical simulations of the dipolar spin ice model [9l [10] and of a dual 
monopole electrolyte. Evidence for the characteristic non-Ohmic conductivity signature 
of a weak monopole electrolyte - the Wien effect - was reported in Refs. |Tll |T2] , where 
the term 'magnetricity' was coined. 

Experimental evidence indicates that magnetic monopoles afford an economical 
description of spin ice at temperatures below ~ 10 K [121 HH US]- In particular, 
down to about 0.3 K the equilibrium specific heat is well described by Debye-Hiickel 
theory [lll|16l|T7]. However, to account in detail for the monopole currents and magnetic 
relaxation is generally a tricky problem, especially in the regime of slow dynamics at 
subkelvin temperatures, and this is an ongoing subject of investigation [121 HH US]- 

Prior to the recent wave of interest, the spin dynamics of spin ice were studied in 
detail by Matsuhira et al. j20] and Snyder et al. [21] by ac-susceptibility, by Ehlers et 
al. [22l|23] using neutron spin echo, by Lago et al. [21] using muon spin relaxation (fiSR), 
by Orendac et al. [26] using bulk magnetocalorimetric methods, by Kitagawa et al. [26] 
using nuclear quadrupole resonance, and by Sutter et al. [27] using nuclear forward 
scattering. Previous studies of the spatial spin correlations by neutron scattering may 
be found, for example, in Refs. [21 EHl [291 |30l EH ESI |33] . 

The spin correlations in the spin ice state are characterised by two remarkable 
features [131 CH EH [35]. The first is a property common to many ice-type models, that 
transverse magnetization (or polarisation) fluctuations are essentially unrestricted while 
longitudinal fluctuations are strongly suppressed at low temperature. This behaviour 
is captured in the phenomenological theory of Youngblood and Axe [36] (formulated 
to describe ice rule ferroelectrics) , in which the deconfined defects do not carry any 
Coulombic charge. The second remarkable property, of course, is that in spin ice 
these defects do carry a magnetic Coulomb charge [H [H]. However, they are also 
associated with a 'Dirac string network' of spin configurations, that while not pairing 
the monopoles, does restrict their motion to some extent [3 [T71 [37] . 

Any complete theory of spin ice needs to account for the difference between 
longitudinal and transverse correlations, the Coulombic interactions of monopoles and 
the effect of the Dirac string network. However, different experiments may pick out one 
or another of these three properties, so approximate models are useful. If monopoles 
are the focus then it is of most interest to discuss the longitudinal response as this is a 
highly sensitive measure of the spatial and temporal correlation of magnetic monopoles, 
as shown below. The simplest approach to treating the monopole correlations is a 
'magnetolyte' model of freely diffusing magnetic charges, in which the effect of Dirac 
strings is subsumed into the transport coefficients, and electrolyte properties such as 
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Debye-Hiickel screening [T7j, Bjerrum pairing [HI UHl [17] and the Wien effect [HI IT2] 
may be naturally formulated. Another (and earlier) approach [6] accounts for the ignored 
spin degrees of freedom in the form of an effective reaction field, and this approach has 
recently been developed to include magnetic charge screening [39] . 

The aim of the present work is to calculate a generalised longitudinal susceptibility 
x(q, cu) for magnetic monopoles in spin ice and to explore its application to experiment. 
The theory described here is only a modest extension of the earlier approaches of 
Ryzhkin [6] and Castelnovo et al. [17], but one that is necessary for the purpose of 
comparing theory with experiment. A useful by-product of this work is a clarification of 
the relationship between these two approaches, and their relation to that of Youngblood 
and Axe [36] . The equations discussed here are valid at temperatures that are sufficiently 
high to avoid the non equilibrium physics of the Wien effect [|] for magnetic monopoles 
(> 0.4 K for Dy2Ti207) [H], but not so high that high energy relaxation processes 
become important (> 10 K for Dy2Ti207) [17]. The results of the present paper are 
applicable to zero and weak applied field {hqH -C 1 T) [§} 

2. Ryzhkin's Approximation 

Ryzhkin [§\ explored the monopole dynamics of spin ice by applying the Jaccard theory 
of water ice defects in the water ice-spin ice analogy. He showed that the magnetic 
current density J = J+ + J_ is related to the rate of change of magnetization M by the 
equation (in our notation): 



and he derived the rate of entropy production associated with the flow of the magnetic 
charges: 



where S is entropy and H is magnetic field. He finally used these relations to derive a 
thermodynamic equation of motion: 



Here k = ucQhq is the monopole conductivity, m = u+ = — u_ is the monopole mobility, 
c = c+ -|- c_ is the total concentration of free monopoles and Q = Q+ = —Q- is the 
monopole charge. The isothermal susceptibility is predicted to be xt = '^Xc where xc 
is the nominal Curie susceptibility for the spin ice system {xc ~ 3.95/T for Dy2Ti207). 

I In spin ice the direct current Wien effect takes the form of a transient increase in charge density 
under apphed field, before equilibration to a value lower than the zero field value |43j : the alternating 
current Wien effect should take the form of a steady state increase in current. 

§ A detailed discussion of the dynamical susceptibility of spin ice in the vicinity of the critical point 
induced by a field of hqH w 1 T applied along the [111] axis ^41, has been given by Shtyk and 
Feigel'man 




irreversible 




(2) 




(3) 
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Eqn. [3] contains much physics and is deserving of subtle appreciation. The term in 
H represents the normal drift current of the charges Q± in the applied field H, and if 
there were only this term, then spin ice would be represented as a true conductor, 
precisely equivalent to an electrolyte. However the term in M opposes the direct 
current and indeed extinguishes it completely at equilibrium, where M = XrH. This 
reaction field does not originate in the magnetic monopoles themselves but rather in the 
configurational entropy of the monopole vacuum: magnetization of the system reduces 
that entropy and hence provides a thermodynamic force that opposes the current [6]. 
It should also be noted that what stops the current in Ryzhkin's formulation is not 
the sample boundaries: this is correct under the approximation that the system is 
homogenous and linear. If one further allows the competition of diffusion and drift to 
set up charge density gradients then the boundaries immediately become relevant and 
one must additionally consider boundary conditions that do not allow the passage of 
monopoles. However this is not necessary in the approximation considered. Finally it 
should be emphasised that the extinction of the current implied by Eqn. |3] holds only 
for very small field and magnetization, for it is only in this limit that Eqn. [2] is valid. 

At first sight the right hand side of Eqn. |3] is zero but this is only true at infinite 
time. By introducing a frequency Fourier transform of the magnetization and field and 
combining Eqn. [T] with Eqn. |3| Ryzhkin found that 



where the inverse relaxation time = k/xt- Eqn. |3]can also be integrated to predict 
the magnetic response to the sudden application or removal of a uniform field (see 
Section. 12). Assuming an ellipsoidal sample, when a uniform external field Hext is 
applied, the internal field Hmt is reduced by the demagnetizing field "DM, such that 
Hint = Hext — ^^M. In spin ice the demagnetizing field arises from the magnetic pole 
density associated with uncompensated surface monopoles. As the spin ice sample is 
magnetized, an imbalance of magnetic charge develops at opposite surfaces as a result 
of the transient monopole current described by Eqn. |3} However, as a result of the 
entropic 'reaction field' discussed above, the monopoles are not sufficiently free to achieve 
complete screening of the internal field. The incomplete screening of of the internal field 
due to magnetic monopoles has been discussed in detail by Ryzhkin and Ryzhkin [39] . 



3. Definition of the two characteristic rates u and uq 

We may define the relaxation rate u by: 

V = T = — = — , (5) 

XT VoXt 

where the concentration c has been substituted for the total dimensionless monopole 
density or mole fraction x = cVq, and Vq is the volume per site of the diamond lattice 
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inhabited by the magnetic monopoles: Vq = (8/3v^)ct^, where a is the near neighbour 
spacing on the diamond lattice. 

Using the Nernst-Einstein relation, the mobility u may be replaced by the diffusivity 

" = It ' 

and hence 

(xt) 

where is the Debye length 

In turn D is determined by the monopole hopping frequency uq. In a simple lattice 
diffusion approximation we may write 



Id. (7) 



D-"^, (9) 

where a is the diamond lattice constant (the numerical factor of 6 may be modified 
very slightly when the fact that a monopole may only hop in three out of four local 
directions is accounted for |.17j). Using the definitions Q = 2/i/a where fi is the rare 
earth magnetic moment, and xc = Ato/^^/3VofcT, Eqns. [sj [6] and |9] may be rearranged 
to give: 

u = guox, (10) 

where (7 = xc/xt = 1/2 in Ryzhkin's theory but is more generally weakly temperature 
dependent and varies between g = 1 at high temperature and g = 2 at low 
temperature |10]. Equation 10 will be seen to be very important for the interpretation 
of experiments on spin ice. 



II The Nernst-Einstein equation, which may be derived from the Boltzmann transport equation, is one 
of the basic equations of electrochemistry. According to Wannier [44], it was used by Nernst in 1888 to 
make the first direct measurement of the elementary electronic charge - at a time when the electron had 
not yet been identified and even the existence of atoms or ions was controversial. The Nernst-Einstein 
equation may be used in this way because diffusivity and mobility are independently measurable for 
an electrolyte. However, for magnetic monopoles in spin ice no way has yet been identified to measure 
these quantities independently. 

^ In the temperature range considered here (e.g. 0.4 - 10.0 K for Dy2Ti207 ), at zero applied 
field, the Debye length is is the dominant length scale in the system, and as shown below, magnetic 
inhomogeneities develop on this length scale. In the lower temperature range, not considered here, 
inhomogeneities on the scale of the Bjerrum length l-r = ^{)Q'^ / ^irkT become particularly important, 
and probably dominate: see Refs. PTIIT^IHB] . 
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4. Coulombic Correlation of the Monopole Current 

We define tlie flux of positive and negative monopoles as j+ and j_ respectively. 
Assuming there is no temperature gradient in the system, then the thermodynamic 
equations of motion are (with i,j = +, — ): 

ji = Lii^i + Lij^j, (11) 

where X denotes a generalised thermodynamic force. If we assume that the monopole 
density is small, then following the theory of weak electrolytes we would expect the 
cross terms Lij with i ^ j to be zero. However the monopole motion should be strongly 
correlated in the sense that it always acts to maintain local charge neutrality: 

Q+j+ + Q_j_ = 0, (12) 

g+c+ + g_c_ = o. (13) 

It is important to emphasise that this thermodynamic force is not the same as Ryzhkin's 
reaction field which is a purely spin phenomenon peculiar to spin ice. In fact, Eqn. |3] 
does not account for Coulombic correlation between magnetic monopoles and in the 
next level of description this needs to be accounted for. 



The left hand side of Eqn. [12] is simply the magnetic diffusion current JdiSusion 
which contributes to the total magnetic current J = SM/St. At equilibrium in zero 
applied magnetic field, the monopole diffusion is such that it does not change the local 
magnetization of the system. Thus positive and negative monopoles tend to move in 
the same direction. After a perturbation, the local magnetic current relaxes to zero, 
even though the monopoles continue to hop around the system. The magnetic diffusion 
current is given by: 

Jdiffusion = E ^di(r) = -DQV6c{v), (14) 

i=± 

where 5c{y) = c+(r) — c_(r). In the zero field equilibrium state the average local gradient 
of charge density is everywhere zero. 



5. Spatial Dependence of Longitudinal Magnetization 

The Coulombic correlations create a diffusion force that tends to smooth the local 
longitudinal magnetization. This may be seen as follows. By Helmholtz' theorem the 
vector field M(r) can be decomposed into an irrotational or longitudinal {i) part and a 
solenoidal or transverse (s) part: 

M(r) = M^(r) + M'(r). (15) 

The spin ice ground state is defined by the condition M*(r) = 0, which gives 

Mground(^) = M"(r). (16) 
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Physically, this is a consequence of the spin ice ground state consisting of closed loops of 
spins. The irrotational or longitudinal part is finite only as a result of thermal excitation 
of magnetic monopoles. As V x M*(r) = 0, the vector Laplacian is simply 

V^M'tr) = V(V-M*(r)), (17) 

a result that will be used below. Henceforth (unless otherwise stated) we shall only deal 
with the longitudinal magnetization and the superscript i will be dropped. 
Defining 0(r) as the magnetic scalar potential, the local magnetic field is: 

H(r) = -V0(r). (18) 

and by Poisson's equation and Maxwell's equation, the local magnetic charge density is: 

QSc{r) = Q(c+(r) - c_(r)) = V ■ H(r) = -VV(r) = -V ■ M(r), (19) 



Thus, using Eqn. [T7| the local charge density gradient is: 

gV5c(r) = -V^M(r). (20) 
The magnetic diffusion current associated with a finite charge density gradient is: 

Jdiffusion = -DQ\/6c{r) = DV'M(r). (21) 

This term can then be added to Eqn. |3] to describe relaxation of the spatial charge 
arrangements: 

J(r) = K (H(r) - Xr^M(r)) + DV^M{r) (22) 

In recent work Ryzhkin and Ryzhkin |39] stated such an equation to facilitate a 
calculation of magnetic screening effects in spin ice. 



6. Free Energy Functional 

The same equation can be derived from a Landau-type free energy functional as follows. 
If we apply a local field H(r) then this induces a longitudinal response M(r). The local 
magnetization is opposed by the entropy cost of ordering the spins of the sample as well 
as the entropy cost of creating a local charge imbalance. Note that these two factors are 
distinct: it is possible to increase order in the sample without creating a local charge 
imbalance. 

From general chemical thermodynamics we expect the entropic cost of charge 
imbalance to make the following contribution to the local Gibbs free energy: 



_fcT (fc(r))^ 
2 c 



local ~ o „ • (2"^) 



Hence using Poisson's equation and Maxwell's equation again: 

GL.= f^^ = f/l,(V.M)^=1^(V.M)= (24) 

/ V Co 2, Zk 



Generalised Longitudinal Susceptibility for Magnetic Monopoles in Spin Ice 



8 



We may then write down a free energy functional for the system: 

G(M(r)) = / ^M{rY - /xoM(r) ■ H(r) + ^ (V ■ M)^ dr, (25) 

where the first term on the right represents the Jaccard entropy [6] which in this 
representation is seen to be equal to the entropy of a cooperative paramagnet |16] . The 
rate of change of longitudinal magnetization M(r) may be found in a linear response 
approximation: 



5M(r) 

dt 



K 



6G(M{r)) 
5M(r) 



H(r) - XT^M{r)] + DV^M{r 



(26) 



which gives Eqn. 22 The derivation of the right hand term of Eqn. 33 is given in the 
footnote from which it can be seen that the contribution of surface charge to the 
Gibbs energy is neglected. 

Under conditions of fixed temperature and field, the rate of dissipation is: 

/dG\ _ _dM{r,t)6G{r,t) 

liicvursiiJic \ / 

Hence, using Eqn. [1} 



'rf) 

V / irreversible 



dt 



6M 



(32) 



irreversible 



D. 



H(r) - Xr^M(r) + — V^M(r 



K 



(33) 



which is the extension of Eqn. [2] to include monopole diffusion. 

Owing to the neglect of surface charge, these equations are generally applicable only 
under conditions of small field and small magnetization, or else at short time. When 
these conditions are violated the build up of surface charge may have a decisive influence 
on the internal fields and on the magnetization process, and the preceding equations 
would need to be modified to account for this. 



+ The functional derivative of G = /(V • M)^c?r is found as follows: 

( "M^, f) - i /{V • (M + ef)rdr = / 2(V • f)(V • M)dr, 



(27) 



where the angular brackets indicate a distribution with vector test function f(r). Now 

V • [(f)(V • M)] = (V • f)(V • M) + (f)(V(V • M)), (28) 
so using Eqn. |17| and the divergence theorem, 
'SG[M{r)] 



(5M 



, f \ = 2 j [(f)(V • M)] -nda- J 2{{){V^M) dr, 



(29) 



where S denotes a surface integral, n a unit normal to the surface and a a surface element. If the 
surface charge is everywhere zero then we have: 

/^G[M(r)] A _ 



and finally 



\ 6M 



6G[M{r)] 



2(f)(V^M) dr, 



(30) 



(31) 
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7. Generalised Longitudinal Susceptibility 

We introduce the spatial and time dependent Fourier transforms of the longitudinal 
magnetization and longitudinal magnetic field: 

M(r, t) = -i- ^ / M(q, io)e^^-^-^-'du, (34) 
ZnV q J 

H(r, t) = ^Y.j H(q, i^)e'''-'-'^'duj, (35) 

and the generalised susceptibility (assuming translational invariance): 

M(q,a;) 



Substituting these definitions into Eqn. 22 and Eqn. [T| we find: 



X(q,a;) = ^. (37) 



where: 



Note that a generalised susceptibility of this sort could also be derived by solving a 
Langevin equation incorporating the Landau free energy, as described in Ref. 



Using Eqn. 10, the generalised susceptibility can also be written: 



X(q,^) = 1 , / 2 2^1 A ■ = c-2 I 2 ^ ■ (39) 

where Tq = I/z/q and the correlation length is: 

It should be emphasised that this is an equation for the longitudinal susceptibility 
only (M II H || q). Whether at equilibrium or not, the latter is finite only if there is a 
finite density of monopoles. In contrast (see Section [sj the transverse susceptibility of 
spin ice could take a paramagnetic value at equilibrium even in the complete absence 
of monopoles, as it does in the monopole-free spin ice ground state. In principle the 
transverse susceptibility could relax through the flipping of closed loops of spins |3], 
though in reality it is more likely that it relaxes through the transient passage of 
magnetic monopoles. A more complete description of the wavevector dependence of 
the susceptibility is given in Refs. [361 HHl El [13], and the possibility of quantum 
mechanical effects giving rise to transverse dynamics distinct from magnetic monopoles 
(so called 'photons') has been discussed in Refs. |l8l H9]. 



There is potentially a problem with Eqn. 38 and the subsequent equations. To see 



this we rewrite i^q as follows: 




'^q = t'o + gx], (41) 
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and assume that the maximum possible equihbrium value of the density is x = 1/2. 
Since it seems implausible that z/q would ever exceed z/q (and indeed uq should generally 
be less than it), then it appears that Eqn. 38 breaks down at g^a^ > 3g. To guarantee 
this never happens we can write: 



1^0 



+ gx 



(42) 



Q{l + Aq^ay3) 

where A is introduced as a phenomenological (dimensionless) parameter. Applying Eqn. 



we find: 



^-2 + g2/(l + Ag2a2/3) - GgiooTo' 



(43) 



As discussed further below (Section 14), it would seem more realistic to use this 
phenomenological equation than Eqn. 39 in order to describe experiment. 



8. Equilibrium Field Fluctuations due to Monopoles 

Consider now the spin ice state in zero applied magnetic field, where the internal 
field H(r), which originates from the magnetic monopoles, is no longer constrained 
but instead relaxes self consistently with the magnetization. An approximation to 
the problem is that close to equilibrium, the field takes the value H(r) = xj^^M(r) 
everywhere, and that this relation is maintained for small fiuctuations away from 
equilibrium. The internal field therefore costs spin entropy: 

-T5field = /io^i^(r)^ (44) 

and energy 

[/field = 2^ ■ H = -Y^'il + Xt), (45) 

but this is offset by the energy gain in magnetizing the sample: 

U' = -fioM ■ H = -/ioXT^'. (46) 

Summing these contributions and allowing the magnetic charge density to fiuctuate 
along with the field we find a free energy functional for field fiuctuations: 

F(H(r)) = / ^H{ry + ^ (V ■ H)^ dr. (47) 

This functional is entirely equivalent to that for electric field fiuctuations in an 
electrolyte and in fact is equivalent to Debye-Hiickel theory |50]. Thus we see that 
by suppressing spin fiuctuations (i.e. setting H(r) = Xj;^M(r)) we recover the Debye- 
Hiickel approximation of Castelnovo et al. [TT] . 

Introducing the Fourier transformed field H(r) = X]q H(q)e*^ '' and 

substituting into Eqn. |47]we find: 

F = ^EHmi + q%\ (48) 
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where V is the volume and we have used D/k = Ij^ and the fact that the field, being 
derived from a scalar potential, is longitudinal to the wavevector. Since the probability 
of a fiuctuation is oc e~^/^^, we immediately see that the mean square amplitude of a 
mode is: 

kT 

( W) = 2,2, ,r (49) 

The mean square field at a point is: 

mrn-uw^-)^j^,i;"j^f^^. (50) 

The integral is dominated by short wavelength modes and the mean square field 
approximately takes the value: 

kT 

{H{rf) ^ . (51) 

9. Mean square induction at a point 



Using Eqn. 51 , the mean square induction at a point is: 



Obviously this equation depends on the induction being averaged over a sufficiently 
large volume that contains locally magnetized spins. If the point with which we are 
concerned experiences no local induction from the magnetized spins, and sees only a far 
field, then the mean square induction is simply 

Recalling that = noQ^x/kTVo where Vq ~ it may be shown that 

{B')i,, ^ Blx (54) 

where 

5o - ^, (55) 

and the symbol ~ is used to indicate that factors of order unity are dropped. In the 
case that there is local induction arising from magnetized spins, this equation may be 
modified to: 

B',-^{1 + Xt). (56) 



It is instructive to derive Eqn. 54 in direct space. At a point in the sample 
the squared field may be averaged over contributions at distance r weighted by the 
probability of finding a monopole at that distance: 

/r2\ ^ /^/^oQ\V4vrr^\ B^x a\ . 

(5)-E(^j ^d^ = Bo^^ (57) 
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which neglects correlation between the field contributions. The fields are correlated over 
a distance of Ijj, but the average number of monopoles within a distance l^ is typically 
of order unity, so correlation may be neglected to a first approximation. The field Bo 
is approximately that due to a monopole or a spin at a distance a. Thus if a defect is 
viewed at a distance a it looks like a spin, but if viewed at a much greater distance it 
looks like a monopole. For this reason monopoles are best detected by measuring their 
far fields [II]. 



10. Relaxation of the Field fluctuations 

The equivalent electrolyte theory has been formulated and worked out in detail by 
Oosawa [50], who found that the relaxation rate of a mode labelled by q is: 

z/q = D{q' + l^') = z/o (^^) + K. (58) 

Thus, short wavelength modes relax at a rate of z/q, the monopole hop rate, and long 
wavelength modes relax at a rate k, the monopole conductivity. Fluctuations are 
important on all scales between the lattice constant and the Debye length, so there 
is a dispersion of relaxation rates from the monopole hop rate uo to the bulk field 
relaxation rate, or monopole conductivity k. 

From Eqn. |50} it may be seen that monopoles at distance a and those at distance 
Ijy make similar contributions to the mean square field, while those at much greater 
distance may be neglected (in zero applied field). Monopoles at distance a reverse the 
local field at a rate of approximately z/q, whereas the cloud of monopoles at distance Id 
only reverses the field at the much slower rate k, although it gives rise to small field 
fiuctuations at a rate uq. 



Comparison of Eqn. 58 with Eqn. |5] suggests that the field correlations relax at a 
rate k while the spin correlations relax at a rate z/ = k/xt- Although this difference 
may reflect the approximations made, it also seems plausible that the spin correlations 
relax more slowly than the fleld correlations at low temperature. Thus the spin system 
can only relax by the passage of monopoles, so the time taken to flnd the most probable 
spin arrangement will generally be longer than the time taken to flnd the most probable 
monopole arrangement. 



11. Spin Lattice Relaxation 

It is obvious that the relaxation rate z/ is equal to the spin-lattice relaxation rate rf^ 
but it is useful to see how this arises in detail. Initial application of a magnetic fleld H 
should result in almost instantaneous magnetization, in which energy is stored within 
the system of effective spins. The spin temperature is therefore initially higher than 
the applied bath temperature T. Following Casimir and du Pre [HH |52] the temperature 
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-ae. 



(59) 



difference 6 = Ts — T determines the rate of exchange of heat Q with the thermal bath, 
and the consequent return of the spin system to thermal equilibrium at temperature T: 
dQ 
dt 

We may also write: 

dQ = CndTs (60) 

where Ch is the heat capacity at constant field. Hence we find: 
dT^ a 



dt Ch 



{Ts - T) 



(61) 



and the spin temperature relaxes at a rate t^^ = a/Cn- 

To link this to the magnetization relaxation we use the thermodynamic relations 

'dH 



dQ = TdS = CndT - fioT 



dT 



xdH, 



M 



and 



xdH = XTdH + 



'DM' 
~dT 



dT, 



(62) 



(63) 



H 



where xdH = dM and Cm is the heat capacity at constant magnetization. If then we 
apply a field H = Hic^^* and elicit a response 9 = Oic^'^^ the above two equations may 
be solved with the substitutions dH — > H,dT — )■ 9, dQ — )■ —aOdT (Eqn. 59) and use 
of the thermodynamic relation Ch — Cm = —fi'oT{dH/dT)M{dM/dT)H- This gives the 
well-known result [52] : 

(Xt - Xs)(^Ti 



X 



Xs + 



XT - Xs 

2^2 



+ 1 



(64) 



where xs = XtCm/Ch is the adiabatic susceptibility. By comparison with Eqn. |4]we 
see that in Ryzhkin's approximation z/ = rf ^ and the adiabatic susceptibility is assumed 
to be zero. The former is easily understood as any magnetization involves the passage 
of a monopole current accompanied by dissipation. However, it is conceivable that the 
adiabatic response could be finite in the real material, and involve the 'stretching' of 
the excited state magnetic moment along the field direction, in which case the adiabatic 
susceptibihty xs could be a direct measure the density of excited states or monopoles. 
This idea needs to be checked in detail. 



11.1. Phonon Bottleneck 

We may also modify this approach to include a 'phonon bottleneck'. The spin system 
is considered to be connected to the phonon system at temperature Tp, and the phonon 
system is connected to the bath at temperature T. The thermal relaxation between 
phonon system and bath is characterised by a thermal conductivity a' . If we make a 
steady state approximation to the phonon temperature, then the rate of heat exchange 
between phonon system and bath is simply: 
dQ dE , , 
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where E = U — fioMH is the magnetic enthalpy. Under the approximation that the 
monopole internal energy U is constant, we find simply that: 

rp = r + (a')-Vo^('^). (66) 

Thus a thermometer placed on the sample could be used to measure Tp and hence gain 
an alternative measure of the magnetic current dM/dt after transients have died away. 



The rise in temperature of the sample (Eqn. 66) occurs when the rate of fiow of 
heat between the spin and phonon systems exceeds the rate of fiow of heat from the 
phonon system to the bath. In the steady state approximation the criterion for this is 
a' < a = uCh- In the low temperature limit we find Ch ~ (|/^|/^r^Vo)a: [16], where 
is the monopole chemical potential [H [7] and hence: 




4^. (67) 

is a criterion for observation of this effect. The ratio x(T, H = 0)/T is always sufficiently 
small that this analysis suggests that the bottleneck can never be observed in zero 
applied field, and hence any observation of a bottleneck is likely to refiect a significant 
field-induced increase in x{T,H) (the Wien effect). This conclusion is consistent with 
the experimental observations of Slobinsky et al. [53] who observed a phonon bottleneck, 
albeit in fields much stronger than those appropriate to the theory discussed here. 

11.2. Thermal Quench 

If bound monopole pairs equilibrate sufficiently quickly with the monopole vacuum, then 
the magnetic monopoles may be regarded as in direct equilibrium with the vacuum: 
(0) = (+) + (— ). The equilibrium constant is 

Keg = Xl (68) 

where we temporarily label the equilibrium density as xq. By definition the 
thermodynamic equilibrium constant is given by K^q = e^'^^'^ where /i < is the 
monopole chemical potential [HIZ]- 

Neglecting Bjerrum pairs [12], the kinetic rate equation for the change in monopole 
density is: 

dx , , , , 

(69) 




dt " V 4 

where the first term on the right accounts for monopole recombination and the second 
for monopole generation. The recombination rate constant has been assumed to be 
equal to the monopole hop rate. If the temperature is lowered at a rate dT/dt = —r 
then it follows that 

^ = -or-^(^-e^'^/^). (70) 

Numerical solution of this equation shows that monopole density reaches a finite 
approximate steady state of the order x = t/uq at low temperatures. Putting in 
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reasonable parameters for Dy2Ti207 spin ice (e.g. uq = 1000 s~^,/i = —4.6 K), a 
rate of cooling of r = 10"^ K s'^ (about 1 K per day) would result in a residual density 
of about X = 10~^ at temperatures lower that 0.25 K (see Fig. [T]). Arrest of the cooling 
at a base temperature significantly less than 0.25 K, where exp(2yu/T) becomes entirely 
negligible, then results in a very slow power law decay of the monopole density according 
to (see Eqn. ItTI) : 



<t) = (71) 

1 +x(0)z/ot/4' ^ ^ 

and even one day of waiting would barely reduce the density by a further power of 10 
(Fig. [2]). Therefore, with any realisable rate of cooling and time of waiting it is not 
possible to completely rid the system of monopoles on experimental time scales. 





0.4 0.6 
Temperature (K) 



Figure 1. Thermal evolution of the monopole density x according to Eqn. 70 (red 



line) versus the equilibrium density 2exp(/i/T) (blue line), with the cooling process 



starting at 1.0 K at a rate 10 ^ K s ^ (here vq 



10^ 



-4.6 K). 



The preceding analysis neglects many factors that may become important at low 
temperatures, including possible thermal evolution of the hop rate, extrinsic factors 
and kinetic constraints arising from the Dirac strings. However, most of these factors 
will tend to tend to reduce, rather than increase, the rate of relaxation, so it is safe 
to conclude that the analysis is correct in its conclusion that a monopole-free state in 
zero applied magnetic field remains inaccessible to experiment A detailed analysis of 

* One may contrast the case of pure water at room temperature, where a density of ^ 10^^ ions per 
water molecule is easily maintained at equiliibrium on experimental time scales: however the diffusion 
constant of in water is some 10^ times larger than that of magnetic monopoles in spin ice. 
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Figure 2. Temporal evolution of the monopole density x starting at x 
according to Eqn. 71 (parameters as in Fig. IT]). 



10- 



idealised thermal quenches in spin ice [37] has identified the important role of monopole- 
antimonopole pairs that cannot immediately annihilate by a single spin flip. These 
'noncontractable' pairs form long lived metastable states at low temperature. 

It is also pointed out in Ref. [37] that at sufficiently low temperatures, owing to the 
divergence of the mobility as 1/T (Eqn, [6]), monopoles will recombine at the maximum 
speed allowed by the monopole hop rate. However, in the electrolyte theory, this effect 
is accounted for by the concept of the Bjerrum pair [TT| [T2]. Such nonlinear response 
occurs only within the pair, that is when the monopole-monopole separation is less 
than It = ^qQ"^ / SnkT . This fast nonlinear response then appears like the flipping of 
giant dipoles of magnetic moment QIt [12], but for monopoles at greater separation. 



the ordinary recombination kinetics of Eqn. 71 are obeyed. The average monopole 
separation grows with decreasing temperature much faster than the Bjerrum pair radius, 
so in a 'slow' quench of the sort described above, the divergence of the mobility should 
not signiflcantly speed up the rate of recombination. The role of Bjerrum pairs, which 
is closely connected to the Wien effect [TTl [T2] is not considered further here. 



12. Application to Experiment: General 

In the following sections I discuss the application of these ideas to different experimental 
measurements. The equations quoted and derived here should be applicable at 
sufficiently small applied field and at temperatures (> 0.4 K for Dy2Ti207) where 
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the Wien effect is absent, so the dimensionless monopole density x depends only on 
temperature. 

It is important to emphasise that the experimental response in all cases depends on 
x{T) = v/ gvQ. There is a general belief that the monopole hop rate z/q is temperature 
independent 0. Assuming this, x{T) can be calculated by numerical simulation [7], 
by Debye-Hiickel theory [T7] , or approximately inferred from the specific heat [16] . For 
Dy2Ti207 the monopole density is roughly constant below 10 K and decreases rapidly 
as the temperature is lowered below 2 K (see Fig. |3|. The corresponding relaxation 
time therefore shows a plateau between 10 K and 2 K, and increases rapidly as the 
temperature is further lowered [7]. The picture of monopoles hopping at a constant rate 
z/q breaks down at temperatures above about 10 K, where an Orbach-like relaxation 
process involving an excited crystal field level becomes important [22]. 

0.6i 1 



0.4 

X 

0.2 
0. 




Temperature (K) 

Figure 3. Approximate evolution of the dimensionless monopole density x with 
temperature [Mj for parameters appropriate to Dy2Ti207. The values of x have been 
inferred from fitting experimental specific heat data to Debye-Hiickel theory, according 
to the method of Ref. [H] . 



13. Magnetization Measurements 

13.1. dc- Magnetization 

To treat a bulk magnetization measurement we can set g = in the above equations. 
In any real sample demagnetizing fields need to be accounted for. If we assume an 
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ellipsoidal sample and write Hint = Hcxt — T)'M. where V is the demagnetizing factor, 
then Ryzhkin's equation becomes: 

«:(He.t-M(x^i + I))). (72) 



dt 

For the case of a steady field this equation may be integrated to find: 

M(t) = + M(0)e-*/^ (73) 

Xt +^ 

so the relaxation of the magnetization is purely exponential. It may be seen that the 
susceptibility xt behaves as an effective demagnetizing field and that the apparent 
susceptibihty is 

_ M _ 1 

which tends towards 1/P as T — )■ 0. 

This equation may be used to describe field cooled (FC) and zero field cooled (ZFC) 
magnetization measurements. It is assumed that in the FC experiment, the sample is 
cooled sufficiently slowly that it always remains in equilibrium (although we have shown 
that this cannot be strictly true), but that in the ZFC experiment it is heated at a 
sufficient rate to be observed on a timescale fobs ^ ^, With these approximations the 
FC and ZFC magnetizations are: 

MzFC = \ (75) 

Xt +^ 



Xt'+T^' 

Using reasonable parameters, these equations predict a large FC-ZFC splitting in 
M{T) / H, as shown in Fig. |4j Here it has been assumed that there is a single observation 
time of about 100 s, which must be a rather crude approximation. Nevertheless, 
a dramatic FC-ZFC splitting, qualitatively similar to that shown, was observed in 
experiment by Snyder et al. [21]. There appear to be two principal ways in which the 
experimental result differs from Fig. |4} First, the experimental FC magnetization below 
the splitting temperature, becomes temperature independent at a value smaller than the 
theoretical M = H/T) [SS]. Second, the experimental splitting temperature (~0.65 K) 
for Dy2Ti207 is higher than that which can be reasonably justified by Ryzhkin's model. 
The higher than expected splitting temperature appears to be related to an anomalous 
slowing down of relaxation seen in ac- magnetization [191 ES], as well as in numerical 
simulations [7|. Possible causes of the experimentally observed slowing down include the 
constraints imposed by of the Dirac string network [Tl|37|, thermal coupling effects [S3] 
and a transition in the monopole density [56]. Also, the Wien effect is important in this 
regime and will play a role in the transient response [121 EH]- 

As regards Ryzhkin's prediction [6] that xt = '^Xc, a recent theoretical study [iO] . 
using parameters appropriate to Ho2Ti207 spin ice, has shown that there is a very 
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Temperature (K) 

Figure 4. Field cooled (FC, blue) versus zero field cooled (ZFC, red) splitting 



according to Eqn. 75 76 based on Ryzhkin's theory of monopole current [B]. The 
observation time has been set at tobs = 100 s, the monopole hop rate at i/q = 10^ 
and the demagnetizing factor at V = 1/3. The monopole density has been roughly 
approximated hy x ~ 2e^^^ j (l + 2e^l^\ with /i = —4.6 K, appropriate to Dy2Ti207. 



slow crossover between xt = Xc at high temperature (~ 100 K) to xt = 2xc in the 
low temperature limit. Experimental measurements appear to be consistent with this 
prediction |10]. This 'Curie law crossover' has not yet been experimentally confirmed 
for Dy2Ti207. 



13.2. ac- Magnetization 

For ac-magnetization measurements, Ryzhkin's equation (Eqn. [s]) can be applied, using 
a demagnetization correction. As described above, the rate u = 1/ri, the spin-lattice 
relaxation rate that arises in the Bloch equations. 

Although Matsuhira et al. have shown that the relaxation is never a simple 
exponential at the temperatures of interest |20j, it appears that the characteristic 
relaxation time does behave according to Ryzhkin's theory. Thus at high temperature 
we would expect a characteristic relaxation time r = l/guox and this is born out in 
experiment in the temperature range > 2 K for Dy2Ti207 where x ~ 1/2 [7]. However at 
lower temperatures, it is evident that the relaxation rate may not be simply proportional 
to the monopole density [3 [191 ES] ■ 
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14. Neutron Scattering 

14- 1- Conventional Neutron Scattering 

Having accounted for the atomic form factor and assuming sufficiently small energy 
transfer, the partial differential cross section of conventional neutron scattering [a") is 
proportional to the imaginary part of the generalised susceptibility: 
kT 

a" oc — Im[x"'^(K.u;)]. (77) 

Here a,P = x,y,z and K is the scattering vector. As shown in Ref. [I3], a 
polarised neutron scattering experiment may be used to isolate the longitudinal (zz) 
susceptibility discussed here by scanning through a Brillouin zone centre perpendicular 
to the reciprocal lattice vector Kq. It is particularly useful to use Kq = (0, 0, 2) in the 
face centred cubic basis as there is no nuclear Bragg peak at that wavevector |13j . 

For scans along this direction (which corresponds to a scan across the "pinch 



point" [T3j, using Eqn. 37 and setting q = K — Kq, we find 



cr' (longitudinal) oc -^—^ ^, (80) 



K,T 

ct" (longitudinal) oc — — ^. (78) 

Unfortunately the dynamics of spin ice are generally too slow to test this expression. 
Instead it is possible to energy integrate and measure in the static approximation 
whereby the differential quasi-elastic cross section a' is given by: 

cr' (longitudinal) oc Tx(q). (79) 

Using Eqn. |37]with u = we find: 

where as already stated, ^ = a/ y/6gx. 

In general x is well approximated by: 
_ exp(-n(T)Jeff/fcT) 

- l + exp(-n(T)JeffAT)' ^ ' 

where Jeg is the effective exchange parameter for a given spin ice Here the prefactor 
n — )■ 4 in the low temperature limit and n — )• 2 in the high temperature limit as a result 
of Debye-Hiickel screening [7]. For Ho2Ti207 spin ice Jes/k ~ 1.8 K, so there is a regime 
at intermediate temperature where ^ ~ exp(1.8/T). In the experiments of Fennell et 
al. [13], the neutron data along wavevectors perpendicular to 002 were fitted to the 
sum of a Lorentzian function and a fiat background. The inverse Lorentzian width was 
indeed found to depend on temperature as predicted here (exp(1.8/T)) although its 
absolute value was much larger than predicted. The fiat background was also found to 



depend on temperature according to Eqn. 81 at high temperatures (with n = 2 and a 
correction for 'double charge' monopoles). 
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14-1-1- Possible Explanation of the Discrepancy There are two potential corrections 
to Eqn. 80 that we did not consider in Ref. [13]. The first stems from the modification 
of Eqn. 39 to give Eqn. 43, as discussed above. Applying this gives: 

cr' (longitudinal) oc Sr, — n n , ^ ■ (82) 

^ ^ ^ + gV(l + Ag2a2/3) ^ ^ 

The second would account for the wavevector dependent misalignment between M(q) 

and q. However, this is a relatively minor correction and is not considered further 

here. Writing a = (v^/4)afcc where Ofcc is the face centred cubic lattice constant and 

q = \/2{2n /afcc)h, we find (for a scan along hhO) 

cr longitudmal oc — ^ . 83 

^ ' (v^7re/a)-2 + /i2/(l + (A7r2/2)/i2) ^ ' 

Using Eqn. |40] this may also be written: 

/n 1- n (3^a;/47r2)TxT f^.. 

cr (longitudmal) oc , , ^, ——, — ; — ^ , , , ^, . (84) 

^ ^ ^ (3^x/47r2) + /i2/(l + (A7r2/2)/i2) ^ ' 

These expressions produce a lineshape and temperature dependence that is very similar 

to that observed in Ref. [13], in that they incorporate both the apparent Lorentzian 

(making it appear anomalously sharp) and the flat background, and they also predict 

the correct temperature dependence in both cases. It would be interesting to compare 

them in detail to the experimental data. 



14-2. Neutron Spin Echo 

Neutron spin echo measures the intermediate scattering function S(K, t) which is 
proportional to the frequency Fourier transform of the right hand side of Eqn. [78] 
Thus we predict 

S{q,t) ~ KTexp(— z/qt), (85) 

with given by Eqn. [5] A test of this expression would require measuring neutron spin 
echo for scattering transverse to the pinch point, as above. Experiments so far [221 123] 
have integrated over larger ranges of q, including transverse fluctuations, and in a 
temperature range where uq. A temperature independent relaxation rate has been 
observed [221 [23], but for Ho2Ti207 this was several order of magnitude faster than that 
derived by ac-susceptibihty on Dy2Ti207. Notwithstanding a possible variation between 
materials it seems likely that the measured relaxation rate is technique dependent, even 
though its temperature dependence is not. This suggests a high frequency component 
to the monopole response that is not contained in the present approximations. 



15. Muon spin relaxation and rotation 

15.1. Longitudinal Field fiSR 

In a /iSR experiment the muon is self trapped by the lattice distortion it creates. In 
a dense magnetic oxide like spin ice it is therefore prone to distort the local magnetic 
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environment that it is aiming to probe. Despite this, the pubhshed results of /xSR 
experiments are reasonably explained by the monopole model. 

Thus a longitudinal field /xSR experiment on Dy2Ti207 was performed by Lago et 
al. [23], who analysed the long time muon depolarisation rate as a measure of the 
field fiuctuation rate. Hence this should have been a measure of u or 1/ri. The 
temperature dependence of the corresponding relaxation time is indeed very close to 
that expected, and it seems very likely that the experiment was observing magnetic 
monopoles. However the magnitude of the relaxation time was an order of magnitude 
smaller than that inferred from ac-magnetization measurements. This would again 
suggest a high frequency component to the monopole response, as noted above. 

15.2. Transverse Field fiSR 

If a muon implants into spin ice at a site of large local field, then transverse field /xSR 
is an uninteresting probe of the spin ice system. Hence we will assume that the muon is 
at a site of zero local field, either within the sample or exterior to the sample, but near 
the surface. While the assumption of zero field sites within the spin ice sample gives a 
highly consistent description of experiment [HI [12], their existence has been contested 
on theoretical grounds [SZ] and the issue has been debated [SHI ESI ED] • 

At sufficiently high temperature (T >10 K) we might expect the TF-/iSR dephasing 
rate A to give a measure of 1/t2, the spin-spin relaxation rate, which may be specified 
by a BPP type [61j expression: T2^- 



where u' is approximately the spin fiip rate and SB^ is the scale of the fluctuations of 
the field component parallel to the applied field. As the latter term tends to dominate, 
we shall only consider this term from now on. 

In the spin ice regime, where x = y/gvo is the dominating parameter of the system, 
the TF-/iSR response is found to have a form that is unfamiliar in the context of /xSR 
on paramagnets. To explain this it is useful to first consider the dimensional analysis of 
the problem. 

15.3. Dimensional Analysis for TF ^SR 

/iSR theory for a simple paramagnet may be formulated in terms of two parameters: 



A and z/'. Here A = 7y {B'^) where the right hand term is the instantaneous root 
mean square field at the muon site, 7 is the muon gyromagnetic ratio, and v' is the 
relaxation rate of this local field. In terms of dimensional analysis we would say that v' 
and A constitute two governing parameters, both with the dimensions of [1 / time\ . The 
quantity of interest in transverse field /xSR is the characteristic rate of muon dephasing. 



1 



1 l''{{5B'f) 




T2 2Ti 
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A. The formal solution to the problem is: 

A = A/ (^^^ , (87) 

where / is an undetermined function. 

In the slow fluctuation limit z/'/A — > and for A to be finite we have / — )■ constant. 
In the fast fluctuation limit A/z/' — )■ and we expect A — )■ 0. The asymptotic form is in 
fact linear in the small parameter, /(1/e) ~ e. The two solutions thus become: 

Asiow ~ A, (88) 
A2 

Afast ~ — , (89) 

formulae that are often used for the analysis of //SR data. 

These formulae may be rationalised by the following heuristic argument. If the 
field sensed by the muon is approximately static on the muon lifetime, then the muons 
precess in phase at a Larmor frequency jBa where Ba is the applied transverse field, 
but accumulate a phase difference A0 = At in time t. If 1/A is equated with the time to 
dephase by order 1 radian then we obtain A ~ A. If , on the other hand, the field jumps 
randomly at rate u, with jump magnitude A, then the phase difference accumulated 
between flips is A0 = A/u' and the phase undergoes a random walk with end to end 



distance vt^A/u) in time t, yielding Eqn. 89 



The case of spin ice is unusual in that there are three, not two, governing parameters. 
The origin of the third governing parameter is in the thermodynamics of the Coulomb 
gas in the grand canonical ensemble where the monopole number N is the sole extensive 
system parameter. We have defined x = N/Nq as a dimensionless monopole density 
(where Nq is the number of diamond lattice sites) and i^o is the temperature independent 
monopole hop rate. As discussed above, the relaxation rate of the local magnetic field 
is u = gxvQ and we may define a scale for the field Aq, that depends only on fixed 
microscopic parameters. 

The formal solution of dimensional analysis can be written: 

Asiow = Ao/f^,-V (90) 



The physical picture we wish to explore is that low temperature (x — )■ 0) corresponds 
to slow fluctuations, and high temperature (x ~ 1) corresponds to fast fluctuations. 
Taking the slow fluctuation limit uq/Aq, u/Aq^O now does not necessarily eliminate 
zz/z/q = gx from the problem. Whether it does so or not depends on the function /. If 
muons detect monopolar fields only (that is, the longitudinal susceptibility), then we 
would expect A to go to zero as a power law in x, for in the absence of monopoles there 
should be no dephasing. In contrast, in the fast fluctuation limit zz/z/q does drop out 



of the problem and we again recover Eqn. 89 The two solutions appropriate to the 
detection of monopolar fields are therefore: 

Asiow-Aof-V, (91) 

\Z/o/ 



Generalised Longitudinal Susceptibility for Magnetic Monopoles in Spin Ice 24 

Afa.t ~ — . (92) 

Thus in the slow fluctuation limit we expect A — z/^, while in the fast fluctuation 
limit, we expect A — )■ constant. The former is an unusual result in the context of /iSR 
and applies to the case where the muons sense only monopolar fields. 



15.4- TF—fjiSR at Low Temperature 

At low temperature the monopole gas is sparse (x ^ 1) and muons that are close 
to monopoles are rare. The muon experiment acts to some extent as a spectroscopy, 
associating different field contributions with different times of observation. Hence to 
use the average field may not be quite correct. The muon signal at long times measures 
only typical muons, which, are far from magnetic monopoles. The typical distance of a 
muon to a monopole is approximately r* /a ^ x~^^'^ and the field sensed by the muon 
is \B\ Bqx^^^. Since this field is random in direction we get the same result for the 
mean square field as above, but with the exponent 4/3 on x instead of 1. In general we 
might expect the apparent mean square field to be given by the equation (S)^ = BqX^, 
with y ~ 1. 

In this limit the Debye length Id is very large and the conductivity k, is very small. 



Although K scales with x, if ?/ < 2 then 7y {B"^) is always larger than it, and the fields 
are quasistatic (here 7 is the muon gyromagnetic ratio). If we approximate the fields 
as completely static on the muon lifetime, then the muons sense a 2;— component of the 
local field that is of the order of the root mean square field. The field sensed by the 
muons is approximately Bq ± J (fi^) and the muons precess coherently at a Larmor 



frequency ^yB^ but are dephased by the spread in local fields. Introducing A = 7^ (-B^) 
the spread of phases accumulated in time t is 

A4> = tA. (93) 

The dephasing time 1/A is equated with the time taken for A0 to become of order one 
radian, with the result [| : 

A'°" ^ = A^ -fBox^/^ (94) 

Hence A is 

.low T 




= 7^0 . (95) 

ji Henley [52] has considered relaxation functions G{t) for nmr for nuclei in sites of zero local field, and 
in the dilute monopole limit has shown that G{t) = exp {^—n{t/Ti)^^ where /3 is a positive exponent 
of order unity and Tj(i = 1,2) is a characteristic timescale for longitudinal (1) and transverse (2) 
relaxation. The results obtained here (for zero local field sites) are consistent with the /iSR relaxation 
function taking this general form. 
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The muon dephasing function depends on the actual field distribution. However it 
is always of the form: 

P = P{Xt). (96) 

This form (with X oc u) was assumed in Ref. [11] and gave a highly consistent description 
of experiment. Although this applies the current ideas in the Wien effect regime, one 
would expect this to be reasonable on the grounds of the dimensional arguments given 
above. Note also that the method of Ref. [11] is insensitive to the precise form of the 
local field distribution. The typical value of A observed in Ref. [TT] was of the order 10^ 
s~^. For DTO spin ice 7-B0 ~ 10* s~^, so a A of 10^ s"'^ corresponds to x = 10~^ if ?/ = 1 
and the monopole field at a typical muon site is about 10~^ T. The temperature at 
which the monopole density is expected to fall to this value is 0.3 K, which is consistent 
with the observations of Ref. [11] . 



15.5. TF—jjiSR at High Temperature 

In the high temperature limit x becomes of order unity so ^ z/q. Thus as we pass 
from low to high temperature, monopoles hopping at a rate vq located near to the 
muon become increasingly important, but as remarked above, these monopoles cannot 
be distinguished from spins, and we return to a model of spin fiipping at rate z/q. In this 
case the ordinary equations of yuSR apply. 

16. Conclusion 

The main conclusion of the present work is that magnetic monopoles in spin ice 
largely determine the longitudinal response of the system. The sole system variable 
for both static and dynamic response is the dimensionless monopole density a;, which is 
determined in a complex way by the four fixed parameters of the problem: a,Q,fi and 
z/q. In contrast, the transverse response does not directly mirror monopole correlations. 

The main theoretical results of this paper are contained in Eqns. [10 , [22, 25, 38, 40 
57} [SH} |M 



Temporal and spatial correlations are linked by x and the Eqns. \10\ and [40] combine 
to establish a dynamic scaling relation: 

z/ ~ (97) 

with ^ = 2, as would be expected for a problem of Brownian motion. It follows (see 



Eqn. 38 ) that there is a dispersion of relaxation rates on all scales from the monopole 



hop rate z/q to the magnetization relaxation rate u. Some evidence has been noted to 



suggest that field fiuctuations relax more quickly than spin fiuctuations (see Eqns. 38 



58 ) but more work is needed to establish this. 

The exponents u and z defined in this way, and the correlation length ^, are not 
conventional quantities as they refiect monopole rather than spin correlations. The 
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spin correlations obey static scaling in the following sense. The correlation function 
g{r), being pseudo dipolar [361 HH] decays as g{r) ~ r~^. Applying the scaling relation 
g(r) r"^'^"^"*'''^ we find i] = 2. As the susceptibility diverges as l/T, the susceptibility 
exponent 7 takes the value 7 = 1. Applying the scaling relation u = 7/(2 — 77) we 
find z/ is infinite, meaning that the spin-spin correlation length remains finite at all 
temperatures. Thus T = marks an unusual critical point with algebraic decay of spin 
correlations, a divergent spin susceptibility, but a non-divergent spin correlation length. 
It is interesting to observe however, that the monopole correlation length does diverge 
at T = 0. 

Free energy functionals for the magnetization and field fluctuations have been 



derived (Eqns. 25, 47) and shown to relate closely to Eqn. 22, previously stated by 



Ryzhkin and Ryzhkin [39l|. In future work it would be interesting to express these 
as functionals of the density (x) and to add further terms to account for energy 
fluxes in the system, as well as Wien dissociation, which both play a role at low 
temperature [m HSl [18] . 

The generalised susceptibility (Eqn. 38) at the level of Ryzhkin's description [6] 
has been derived, as well as the field fluctuation at the level of Debye-Hiickel theory [T7] . 
The latter was used to calculate the longitudinal field fluctuation at a point in the system 
(Eqns. 53, 54), which may be compared to established results for electrolytes [50]. The 
expressions for the mean square field distribution have been used to show that a point 
probe such as a muon at a 'spin free' site (either inside the sample or just outside) will 
give a direct measure of the monopole density as assumed in Ref. [H]. 

It has been shown that according to the electrolyte theory, a non-equilibrium 
population of monopoles is always frozen into the sample, regardless of the rate of 
cooling. However in sufficiently weak magnetic field there is never a phonon bottleneck. 
The former effect should generally be considered when treating low temperature 
experimental data. 

In general the theory discussed here works qualitatively well for real spin ice 
materials, capturing the temperature, wavevector and time dependence of a diverse 
range of experimental responses. However, there are three clear discrepancies. First, 
while the temperature and wavevector dependence of the neutron scattering cross section 
are well accounted for, the amplitude of the correlation length is not, being an order of 
magnitude longer in experiment compared to theory [13]; a possible explanation of this 
has been proposed here. Second, while different experiments [211 EDI ESI |23l [21] agree 
on the temperature dependence of the relaxation rate, they exhibit a wide range of 
relaxation rates: it appears that there is a high frequency response, not accounted for in 
the hopping model. Third, the ac-susceptibility relaxation in the high temperature limit 
is more strongly dispersed |2i[Ij than predicted by the simple approximations discussed 
here. 

It seems very unlikely that the monopole theory will have to be abandoned to 
explain these discrepancies. More likely it needs to be refined. In addition to the 
possible revision of the neutron scattering line shape discussed above, one might also 
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need to consider the effect of quantum ffuctuations |1HI US] or minor terms in tlie spin 
ice Hamiltonian |32]. Also, microscopic factors affecting tfie rate of local spin flipping 
or monopole hopping probably remain to be identified. However, it should also be 
emphasised that the most distinctive aspect of Coulombic correlation - the tendency 
to Bjerrum pairing - has not been accounted for here, or in other 'high temperature' 
theories, but will certainly play a role. Thus Bjerrum pairs have been argued to be 
important in the low temperature non equilibrium regime [12] and have been identified 
in specific heat measurements p5]. Finally, the Wien effect [63], though weak at the 
'high' temperatures considered here, still exists in a screened form [61], and should be 
accounted for in a more accurate description. Although there is much work to be done, 
it is clear that the monopole theory of spin ice [H [6] is a remarkably simple and effective 
description of a complex condensed matter system. 
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